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Abstract 

We study the fractional gravity for spacetimes with non-integer 
dimensions. Our constructions are based on a geometric formalism 
with the fractional Caputo derivative and integral calculus adapted 
to nonolonomic distributions. This allows us to define a fractional 
spacetime geometry with fundamental geometric/physical objects and 
a generalized tensor calculus all being similar to respective integer di- 
mension constructions. Such models of fractional gravity mimic the 
Einstein gravity theory and various Lagrange-Finsler and Hamilton- 
Cartan generalizations in nonholonomic variables. The approach sug- 
gests a number of new implications for gravity and matter field the- 
ories with singular, stochastic, kinetic, fractal, memory etc processes. 
We prove that the fractional gravitational field equations can be in- 
tegrated in very general forms following the anholonomic deformation 
method for constructing exact solutions. Finally, we study some exam- 
ples of fractional black hole solutions, fractional ellipsoid gravitational 
configurations and imbedding of such objects in fractional solitonic 
backgrounds. 
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1 Introduction 



Modern classical and quantum gravity models are constructed geomet- 
rically to be extremely complex with black holes, singularities and horizons, 
nontrivial topology, stochastic processes and fractals, irregular sets, possi- 
ble extra dimensions etc. Such configurations are defined as exact and/or 
approximate solutions of Einstein equations in general relativity and modi- 
fications. 

There is a recent interest in fractional analysis and non-integer dimen- 
sion geometries which resulted in an increasing number of publications dur- 
ing last decade. Calculus with derivatives and integrals of fractional order is 
applied in various directions in physics, for instance, in Lagrange and Hamil- 
ton mechanics and dynamical models |331 [25j |4"T1 [3] , fractal and chaotic 
dynamics [ITJ [HI [381 EU [87] , quantum mechanics [HI [171 [26] , kinetic the- 
ories [861 [881 S3] i plasma physics [101 [39] and astrophysics and cosmology 
[4Q| 127} [28l [35] , and many additional applications in physics and other sci- 
ences [851 El 13(37]. 

The cornerstone questions to be solved in constructing classical and 
quantum gravity theories are positively connected to the meaning of space- 
time dimension, extra dimension and/or fractional dimensions physical ef- 
fects, possible contributions from (non) holonomic and/or commutative/ 
noncommutative variables, of fractal dimension, nonlocal field theories etc. 
With respect to dark energy and dark matter models, various attempts to 
construct the quantum gravity theory, to determine the status of singular- 
ities in fundamental theories etc, it is thus natural to pursue alternative 
concepts of differential and integral calculus and consider new ideas and 
models of spacetime geometry. 

The fractional dimension paradigm is at present extended to new ideas 
on gauge invariance in fractional field theories |13j . being constructed so- 
lutions of fractional Dirac equation [32 \ 124]. with corresponding formalisms 
for generalized Clifford algebra [89] . A formulation of Noether's theorem 
for fractional classical fields |23] and, recently, a model of fractional gravity 
based on Riemann-Liouville concept of non-integer dimensional derivative 
was proposed in [22J. 

Let us provide a simple motivation why knew mathematical tools from 
non-integer calculus may be very useful in solving fundamental theories in 
physics. An exact solution of Einstein equations with singularity for r — > ro 
(for instance, inducing a singular term in the Riemann curvature tensor) 
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modelled by a function 

/(r) ~ V(ro)(r - r )" * s + 2 f(r )(r - r ) * s + .., (1) 

for some real values /(rn) and 2 /(ro) and positive integers x s and 2 s. 
The singular properties of such a function may have a complete different 
geometric and classical/quantum physical meaning if instead of the standard 
partial derivative d/dr we work with the left Caputo fractional derivative 

a 

ir d r f(r) (when usual partial derivatives are generalized through a more 
complex integral-derivative relation, see below formula ([2]))Q Such functions 
positively result in various new physical implications if there are developed 
models with inhomogeneous, fractal, stochastic, kinetic etc processes, in 
classical and quantum theories. 

The advantages of fractional derivatives and integrals become already 
apparent in modern mechanics and for describing thermodynamical and 
electrical properties of real complex media (as we emphasized in the above 
references). We suppose that for non-integer dimensions, new conceptual 
features of gravitational nonlinear interactions can be analyzed. Fractional- 
order models are more appropriate than the integer-order ones because 
they provide an excellent instrument for describing processes with memory, 
branching and hereditary properties which in the classical (pseudo) Rie- 
mannian approach is, in fact, neglected (even such solutions can be encoded 
into the structure of certain quantum/stochastic modifications of Einstein 
equations) . 

In our recent work [75], we provided some theoretical arguments that 
fractional spacetime theories should receive more attention in modern geo- 
metric analysis, evolution theories and gravity. Such a conclusion was de- 
rived from the facts that the Ricci flow theory can be naturally generalized 
for non-integer dimensions, with a new type of statistical geometric analogy, 
if a corresponding nonholonomic fractional differential-integral calculus is 
applied. Such a scenario may transform (pseudo) Riemannian metrics into 
certain fractional ones, or inversely. This involve a number of consequences 
in classical and quantum theories of fundamental field interactions and evo- 
lution models. Here we also remember that the theory of fractional calculus 
with derivatives and integrals of non-integer order goes back to Leibniz, 
Louville, Grunwald, Letnikov and Riemann [301 l3~T] \15\ [29] (The first exam- 
ple of derivative of order a = 1/2 was considered by Leibnitz in 1695, see 
historical remarks in |36j). 

1 We follow the system of notations with left abstract labels explained in Refs. |75II61| . 
see also next section. 
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One might hope that extending a fundamental physical theory like gen- 
eral gravity to fractional dimensions will solve the most troublesome prob- 
lems of Planckian physics, quantum non-renormalizable theories, new types 
of symmetries etc. Such extensions should be constructive ones allowing ex- 
act solutions and possible simple geometric and physical interpretations. In 
this paper, we shall apply the fractional geometric formalism developed in 
Ref. [75] to formulate and study a version of fractional Einstein type gravity 
theory. Certain constructions with the fractional Ricci tensor, and general- 
izations to osculator bundles (higher order fractional mechanical models etc) 
and non-integer dimension Einstein tensor were provided in Refs. [HE2]. It 
is a very cumbersome task, both technically and theoretically, to solve field 
and evolution fractional equations in exact form following those approaches 
with the Riemann-Liouville (RL) fractional derivative. 

In general, there are various possibilities to derive fractional differential- 
integral structures and the question on which derivative would be more 
suitable for formulating engineering and scientific problems remain an open 
issue. Here we note that integrals on net fractals can be approximated by the 
left-sided RL fractional integrals of functions, see |34| . Such a fractal and 
non-integer calculus formalism is applied for certain attempts to construct 
the quantum field theory, gravity and cosmology in fractal universes and 
related diffusion processes in string theory, see J8[ [9] and references therein. 
Nevertheless, we shall follow an approach based on Caputo fractional deriva- 
tives allowing, for instance, to get nonhomogeneous cosmological solutions 
generalizing the off-diagonal metrics from [74]. 

This paper has three purposes: The first one is to formulate a fractional 
generalization of the Einstein gravitational field equations with nonholo- 
nomic local frames and deformations of connections induced by the Caputo 
left derivative. The second one is to prove that such nonlinear systems of 
partial derivative and integral equations can be solved in very general form 
using our anholonomic deformation method of constructing exact solutions 
in gravity and Ricci flow theories, see reviews of results and applications in 
[531 E3 E3 ED EH [73]. The third purpose, is to study certain applications 
of fractional calculus in black hole physics: we construct in explicit form 
such fractional black hole solutions and study their imbedding in fractional 
solitonic backgrounds. 

The article is structured as follows. In section [2] we provide an introduc- 
tion into fractional Caputo type calculus on nonholonomic manifoldsil The 

2 In "integer" differential geometry, there are used also equivalent terms like anholo- 
nomic and non-integrable manifolds. A nonholonomic manifold of integer dimension is 
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Einstein equations are generalized on fractional nonholonomic manifold in 
section [3l Then, in section H] we show that in our approach the fractional 
field equations for gravity can be integrated in very general forms. We 
provide some applications in modern gravity by considering black hole solu- 
tions in factional spacetimes, see section Discussion of results and future 
perspectives are presented in concluding section [6l 

2 Fractional Derivatives, Integrals and Forms 

We provide an introduction to fractional differential and integral calculus 
on "flat" spaces following [32] and then extend the approach as fractional 
generalizations of (nonholonomic) Einstein spacetimes. Relevant details and 
references can be found also in a partner paper [73] (see there and Appendix 
containing a summary of non-integer calculus and methods). 

2.1 Fractional Caputo derivatives 

There are different approaches to fractional derivative and integral calcu- 
lus when there are used different types of fractional derivatives. We elaborate 
and follow a method with nonholonomic distributions when the geometric 
constructions are most closed to "integer" calculus. 

Let us consider that f(x) is a derivable function f : [ \x, 2%] for 
R 3a > 0, and denote the derivative on x as d x = d/dx. The fractional 
Caputo derivatives are determined respectively by 

left: x J x /(x) := f ^-y J (x - x y~«-i (JLY f( x > )dx >. : (2) 

\x 

right: x d 2X f{x) := | 1 (-^J f(x')dx' , 

X 

where we underline the partial derivative symbol, d, in order to distinguish 
the Caputo operators from the RL ones with usual <9|1 A very important 

a 

property is that for a constant C, for instance, lX d x C = 0. 

defined by a pair (V, TV), where V is a manifold and N is a nonintegrable distribution on 
V, see the next section for a generalization to non-integer dimensions. 
3 The left Riemann-Liouville (RL) derivative is 

:= (sL) / {x - -') s " a ~ i /(- , )d-', 
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2.1.1 Fractional integral 

We denote by L z ( \x, 2%) the set of those Lesbegue measurable functions 

f on [ ix, 2 X \ for which \\f\\ z = (/ \f(x)\ z dx) 1 l z < oo and write C z [ \x, 2 X \ 

for a space of functions which are z times continuously differentiate on this 
interval. 

For any real-valued function f(x) defined on a closed interval [ ±x, 2%], 

a 

there is a function F(x) = lX I x f(x) defined by the fractional Riemann- 

a x 

Liouville integral lX I x f(x) := J (x — x') a f(x')dx', when the function 

ix 

a 1-1 

f{x) = lX d x F(x), for all x £ [ ±x, 2^]J1I satisfies the conditions 

a / a \ 

ixOxi ixlxf(x)j = f(x), a>0, 

ix% ( 1 xd x F{x)^j = F(x) - F( lX ), < a < 1. 

A fractional volume integral is a triple fractional integral within a region 
X C M 3 , for instance, of a scalar field f(x k ), 

1(f) = I[x k ]f(x k ) = V] V] I[x 3 ]f(x h ). 

a 

For a = 1 and f(x,y,z), 1(f) = fff dVf(x,y,z) = fff dxdydz f(x,y,z). 

x x 

2.1.2 Fractional differential forms 

An exterior fractional differential can be defined through the fractional 
Caputo derivatives which is self-consistent with the definition of the frac- 
tional integral considered above. We write the fractional absolute differential 

where F is the Euler's gamma function. The left fractional Liouville derivative of order a, 

ct a 

where s — 1 < a < s, with respect to coordinate x is d x f(x) := lmijaj-^-oo iX dxf(x). 
The right RL derivative is x d 2X f(x) := (~^Y j (x' - x) 3 -"- 1 f(x')dx' . The 

X 

a ol 

right fractional Liouville derivative is x df(x k ) := lim 2a; _ >00 x d 2X f(x). 

ct 

The fractional RL derivative of a constant C is not zero but, for instance, lx d x C = 

C — • Integro-differential constructions based only on such derivatives seem to be 

very cumbersome and has a number of properties which are very different from similar 
ones for the integer calculus. 

4 this follows from the fundamental theorem of fractional calculus [42] 
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a 

d in the form 



d := (dx 3 ) a ofy, where dx 1 = {dx ] ) a 1 ' 



r(2-a)' 

where we consider \x % = oil 

An exterior fractional differential is defined 

n 

d = r(2 - a)(x J ') a-1 dar* oij- 
i=i 

Differentials are dual to partial derivatives, and derivation is inverse to inte- 
gration. For a fractional calculus, the concept of "dual" and "inverse" have 
a more sophisticate relation to "integration" and, in result, there is a more 
complex relation between forms and vectors. 

The fractional integration for differential forms on an interval 
L = [ \x, 2x] is defined 

a a 

L I[x] lX d x f{x) = /( 2 x) - /( ix), (3) 

a 

i.e. the fractional differential of a function f(x) is lX d x f(x) = [...], when 

The exact fractional differential 0-form is a fractional differential of the 
function 

a a 

lX d x f(x) := {dx) a lX d x ,f(x'), 

when the equation ([3]) is considered as the fractional generalization of the 
integral for a differential 1-form. So, the formula for the fractional exterior 
derivative can be written 

lX d x := (dx l ) a ix d { . (4) 

For instance, for the fractional differential 1-form oo with coefficients {uji(x k )} 
is 

£ = {dx l ) a iOi(x k ) (5) 



For the "integer" calculus, we use as local coordinate co-bases/-frames the differentials 
dx j = (dx j ) a=1 . For < a < 1 we have dx = (dx) 1 ~ a (dx) a . The "fractional" symbol 

(dx J ) a , related to dx 3 , is used instead of dx' for elaborating a co-vector/differential form 
calculus, see below the formula ((6]). 



8 



and the exterior fractional derivatives of such a fractional 1-form ui gives 
a fractional 2-form, lX d x ( u) = {dx l ) a A (dx 1 ) 01 ix dj Ui{x ). Such a rule 

a 

[15] follows from the property that for any type fractional derivative d x , 
we have 

ot ^ „ / rv \ f a—k „ \ a=k 

dA 1 / 2 /) = E( k )( 9 * f) d * ( V), 

fc=o V / \ / 

when for integer fc, 



a 



7 and <9 X = 0, ft; > 1. 



fc / r(l - a)T(k + 1) 
In fractional derivative calculus, there are used also the properties: 



, x 'd x / ix' — \x'Y = t^tt, KH 



r(/3 + i) 



CS + 1 - a) 



[x — ix) 



where n — 1 < a < n and (3 > n, and ix >d x / [x 1 — \x') = for ft; 



0, 1, 2, n — 1. We obtain 



^(X - XX)" = (dx)" = (dx)°r(Q + 1), 

i.e. 

^ a = r( \u ^ x ~ ix)a - (6) 
T(a + 1) 

We might write the fractional exterior derivative ([4]) in the form 



1 

T(a + 1) 



\xd x '■ — 7t \ X d x (x \x ) 



and the fractional differential 1-form ([5]) as 

a 1 a . . 

T(a + 1) 

Having a well defined exterior calculus of fractional differential forms 
on flat spaces M. n , we can generalize the constructions for a real manifold 
M, dim M = n. 
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2.2 Fractional manifolds and tangent bundles 

A real manifold M, with integer dimension dimM = n, can be endowed 
on charts of a covering atlas with a fractional derivative-integral structure 
of Caputo type as we explained above. In brief, such a space (of necessary 

a 

smooth class) M will be called a fractional manifold. 

Let us explain how the concept of tangent bundle can be developed for 
fractional dimensions. A tangent bundle TM over a manifold M of integer 
dimension is canonically defined by its local integer differential structure <9j. 
A fractional generalization can be obtained if instead of d% we consider the 

left Caputo derivatives ix id_i of type (|2j), for every local coordinate x % on 

a 

a local cart X on M. A fractional tangent bundle TM for a G (0, 1) (the 
symbol T is underlined in order to emphasize that we shall associate the 
approach to a fractional Caputo derivative), see more details in |75|FI For 
simplicity, we shall write both for integer and fractional tangent bundles the 
local coordinates in the form u 13 = {x? ,yi). 

a 

On TM, an arbitrary fractional frame basis 

h = zV)%> (7) 

is connected via a vierlbein transform e /3 g (u' 3 ) with a fractional local coor- 
dinate basis 

a /a a a a \ 

dw = d f = , d f ,d b , = y d b , (8) 



& - XHj' - lX J> Vjl,V b > ~ iy b 

when j' = 1, 2, n and b' = n+ 1, n + 2, n + n. There are also fractional 
co-bases which are dual to (J7]), 

% P = e${vPjdvP, (9) 
where the fractional local coordinate co-basis 

= ((dx i 'r,(dy a 'r), (io) 

when the h- and v-components, (dx l ') a and {dy a ') a are of type ([5]). For 

of o 

integer values, a matrix er a is inverse to eJj, ■ 

a a 

Similarly to TM, we can define a fractional vector bundle E on M, when 
the fiber indices of bases run values a' , b' , ... = n + 1, n + 2, n + m. 



We to not follow alternative constructions with RL fractional derivative [TJ[22] because 
in that direction it is not clear how to prove integrability of the RL-fractional Einstein 
equations. 
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3 Einstein Equations on Fractional Manifolds 

We provide a fractional generalization of the Einstein gravity of arbitrary 
dimensions via nonholonomic variables in such a form when we shall be able 
to construct exact solutions of the gravitational field equations. 

3.1 Introduction to the geometry of fractional manifolds 

Let us consider a "prime" (pseudo) Riemannian manifold V is of inte- 
ger dimension dim V = n + m,n > 2,m > 1. Its fractional extension is 

a 

modelled as a fractional nonholonomic manifold V defined by a quadruple 

a a a a 

(V, AT, d, I), where TV (see below formula (jlip ) is a nonholonomic distri- 
bution defining a nonlinear connection structure, the fractional differential 

a a 

structure d is given by $7$) and Q and the non-integer integral structure I 
is computed following the by rules of type Q. 

3.1.1 Nonholonomic fractional distributions and frames 

a a 

A nonlinear connection (N-connection) N for V is defined by a nonholo- 
nomic distribution (Whitney sum) with conventional h- and v-subspaces, 

a a 

hV and vV, 

a a a a 

TV = hV®vV. (11) 

a 

Nonholonomic manifolds with a N are called, in brief, N-anholonomic 
fractional manifolds. Locally, a fractional N-connection is defined by its 

a 

coefficients, N={ a iV~"}, when 



a 



N= a N?(u)(dx l ) a ®da- (12) 

a 

For a N-connection N, we can always a class of fractional (co) frames (N- 
adapted) linearly depending on a Nf), 



3 rf 



(13) 



a e p = [* e J = (d x 3) a , a e b = (dy b ) a + a N%(dx k ) a ]. (14) 

a 

The nontrivial nonholonomy coefficients are computed a Wf h = d b a N? and 
a W§ = = a e { a Nf - a ej a Nf (where a fl^ are the coefficients of 

the N-connection curvature) for 

[ a e a , a e p ] = a e a a e p - a e p a e a = a W^ a e 7 . 
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3.1.2 Fractional metrics and distinguished connections 

a a 

A (fractional) metric structure g = { a g a /3} is defined on a V by a 
symmetric second rank tensor 

g= a g2l (u)(duir®(duP) a . (15) 

For N-adapted constructions, it is important to use the property that any 

Ct 

fractional metric g can be represented equivalently as a distinguished metric 
(d-metric), g = [ a g kj , a g cb ] , when 

g = a g kj (x,y) a e k ® V + a g cb (x,y) a e c a e b (16) 
= Vk'f °e k ' ® a e j ' + V Cb> °e c ' ® V , (17) 

where matrices r^y = cfzagfil, ±1, ±1] and r/ a /f,/ = diag[±l, ±1, ±1] 
(reflecting signature of "prime" spacetime V) are obtained by frame trans- 
forms 

Vk'f = e k k> e J j> a 9kj and r] a >v = e a a , e b b , a g ab . (18) 

For fractional computations, it is convenient to work with constants ij^/ji and 
r\ a i b i because the Caputo derivatives of constants are zero. This allows us to 
keep the same tensor rules as for the integer dimensions even the rules for 
taking local derivatives became more sophisticate because of N-coefficients 
a Nf(u) and additional vierbein transforms e k k ,(u) and e a a , (u). Such coef- 

ficients mix fractional derivatives computed as a local integration ([2]). 
If we work with RL fractional derivatives, the computation become very 
sophisticate with nonlinear mixing of integration, partial derivatives etc. 

a a 

A distinguished connection (d-connection) D on V is a linear connec- 
tion preserving under parallel transports the Whitney sum (jlip . Using the 
formalism of fractional differential forms introduced in previous section, we 
can elaborate a covariant fractional N-adapted calculus on nonholonomic 

a 

manifolds. To a fractional d-connection D we can associate a N-adapted 
differential 1-form of type ([5]) 

a T T p = a T T Pl a e\ (19) 

when the coefficients are computed with respect to (|14h and (|13p and para- 
metrized the form a T\ p = ( a L) k , a L a bk , Q Cj c , a C^ . 

a 

On fractional forms V, we can act with the absolute fractional differential 

a a 

a d = lX d x + xydy. In N-adapted fractional form, the value Q d := a er "e^ 
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consists from exterior h- and v-derivatives of type , when 

lX d x := {dxT 1 X = a eP a ej and iy d y := {dy a ) a 1 J a = a e b a e b . 
3.1.3 Torsion and curvature of fractional d connections 

a 

The torsion and curvature of a fractional d-connection D = { a T T p } 
are computed, respectively, as fractional 2-forms, 



a T T = D a e T = a d a e T + a T T A a e? and (20) 



ft/3 D r /3 = d T^- T'^A r 7 = R^ e r A e. 
The fractional Ricci tensor a lZic = { a R Q/ g = a R T a/3T } is 

a r> — a 7?^ Q F? — _ a a f? — a F?k Q /? — a Tf c f91 "\ 

a 

The scalar curvature of a fractional d-connection D is 

<*R= « g T^« Rr ^ = QjR+ a S) <* R = a g ij aR .^ a S= a g ab a R ^ 



defined by a sum the h- and v-components of (12ip and contractions with 
the inverse coefficients to a d-metric (|16p . 

We can introduce the Einstein tensor a £ns = { a G a p}, 

a G a/ 3 := "R a /3 — - a g a /3 s (22) 

This allows us to elaborate various types of fractional models of gravity (for 
different types of d-connections and fractional matter sources) generalizing 
the Einstein gravity theory and various modifications. 

Finally, we note that for integer values of a the above formulas transform 
into similar ones for nonholonomic manifolds [57\ \61\ [70] . 

3.1.4 Preferred fractional linear connections 

For applications in modern geometry and standard models of physics, 
there are considered more special classes of d-connections: 

a 

• On a fractional nonholonomic V, there is a unique canonical fractional 
d-connection a D = { a f\ p = ( a L) k , a L a bk , a C) c , a C^} which is 
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compatible with the metric structure, a D ( a g) = 0, and satisfies the 
% fc = 0and «f% c 



conditions = and a T a kr , = 0@ 



• The Levi-Civita connection a V = { °T 7 a ^} can be defined in stan- 
dard from but by using the fractional Caputo left derivatives acting 
the coefficients of a fractional metric (1151). 



As a consequence of nonholonomic structure, it is preferred to work on 

V with a D = { a T T a} instead of Q V (the last one is not adapted to the 

N-connection splitting (jlip ). The torsion a T T (|20p of Q D is uniquely 
induced nonholonomically by off-diagonal coefficients of the d-metric (116p . 

With respect to N-adapted fractional bases (fl~3)) and ([HP , the coeffi- 
cients of the fractional Levi-Civita and canonical d-connection satisfy the 
distorting relation^] 

a ?\p= a T\ P + "^V (23) 

It is not possible to get relations of type ([23]) if the fractional integro- 
differential structure would be not elaborated in N-adapted form for the 



7 The N-adapted coefficients are explicitly determined by the (|16[1 . 

a*ri 1 a ir / a a . a a a a \ 

L jh = g 9 ( Gk 9jr + G] 9hr ~ Gr 9jk >' 

ct "r"a ot / ct \r<i\ i ot ac ( at Ot ot ct ot \rd ot a ct \rd\ 

L bk = e b ( N k j + - g e k g bc - g dc e b N k - g db e c N k J , 

ot s~-yi 1 ot ik ot ot ot s~ia 1 ot ad / ot ot , ot ct ot ot \ 

Cjc = 2 9 6c * fe ' bc= 2 9 ' 6c SM 6c 5cd ~ 6d 9bc ' ' 

We can verify that introducing above formulas into (|20[1 we obtain that T l j k = and 
T^, c = 0, but T l j a , T a jj and T a bi are not zero, and that the metricity conditions are satisfied 
in component form. 

8 the N-adapted coefficients of distortion tensor Z 1 B axe computed 

ct rya ct ct a ab 1 ' ' ( )" 

^jk — — <^j b 9ik 9 — 2 i fc ' 

ct r^c a ct ii 1 / si ch ct ct ih\ ct/~tj 

Mjk 9cb g J -^{SjS h - g jk g ) C J hb , 

t ra cb i ot ct ab\ r ot r c ot / ot atcni 

(d c d d + g cd g )[ L bk - e b ( N k )\ , 

ct r^a ot a ii ■ 1 / ei ch ot ot ih\ a s~lj 

U jk g C b g J +^{0jS k - g jk g ) C J hb , 

/ ra cd a ot ad\ [ ot r o ot / a Arc\l a rrd r\ 

-^(o c o b - g cb 9 )[ L dj - e d { Nj)\, Z bc = 0, 

a ij 

y f [ ol t c a / a »tc\1 ct i V ct t c ct / a atc\1 ct t 

— — {[ L aj - e a ( N 3 )\ g cb +[ L bj - e b ( N 3 )\ g ca \. 



Ot ryi 


= o, 


Ot ryi 

^bk 


1 


_ 2 


Ot iya 

^bk 


1 

_ 2 


Ot ryi 

^kb 


1 

~ 2 


Ot rjCL 

*jb 




a ryi 
^ab 
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left Caputo derivative. 

3.2 Fractional Einstein equations for connections D and V 

In this section, we show that the fractional gravitational equations with 
Caputo fractional derivatives can be integrated in general form similarly to 
the results for integer dimensions [5U EH E3 ED ED] . 

3.2.1 Nonholonomic variables in general relativity 

The Einstein equations on a spacetime manifold of integer dimension V, 
for an energy-momentum source of matter T a g, are written in the form 

E/3-y = R ps- T^gpsR = xTps, (24) 

where k = const and the Einstein tensor is computed for the Levi-Civita 
connection V. It is not possible to integrate in any general form such non- 
linear systems of partial differential equations. 

The Einstein equations (|24p can be rewritten equivalently using the 
canonical d-connection D = {r 7 Q o}, 

Eps = Rps - ~gps S R = Yps, (25) 
L c aj = e a (N?), C} b = 0, n a ji = 0, (26) 

where R p$ is the Ricci tensor for r 7 Q/3 , S R = g^R ps and is such 
a way constructed that Tps >cTp$ for D — >■ V. In general, the Einstein 
tensor E p$ in ([25]) is not equal to the Einstein tensor Ep 1 in (p4|). 

There are two possibilities to make equivalent two different systems of 
equations for V and, respectively, for D. In the first case, we can include 
the contributions of distortion tensor Z^ a a from (|23[) into the source T ps ~ 
xTps+ Z Y pslZ^p], in such a form that the system (|25l) is equivalent to (|24l) 
(both types of such systems of equations are for the same metric structure 
g^ but in terms of different connections). In the second case, we consider 
that Yps = xTps but in order to keep fundamental the Einstein equations 
for V (even for some purposes we shall prefer to work with D) we have to 
impose (at some final steps) the constraints (I26D when the tensors TT^ ([20l) 

and Z 1 a /; (j23"j) are zero. For such constraints, we have r 7 a/3 = T 7 a/3 , with 
respect to N-adapted frames (fl~3|) and (fl~4"|) . even D ^ V. 
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It is very surprising that for some general ansatz for metrics (see below 
the integer version of (|29|) ) there is a separation of equations (|25p for D which 
allows us to integrate such systems in very general forms. If we work from 
the very beginning with V, it is not possible to get from (|24p a generally 
solvable system of equations. Our idea was to encode the geometric and 
physical data for V into D and work with N-adapted constructions when, 
for instance, certain general solutions can be found "easily" for various types 
of non-Riemannian and/or Lagrange-Finsler theories. Imposing additional 
constraints (I26D . we are able to extract Levi-Civita configurations with V, 
for instance, in Einstein gravity. This way, see original results and reviews 
in a series of our works [2U \57\ EU EH [70] , we elaborated a new method 
of constructing exact solutions with generic off-diagonal metrics in gravity 
theories (the so-called anholonomic deformation/frame method). Such a 
geometric technique seems also to be very efficient in Ricci flow theories |62l 
[63| l65| l66| 167] l68| [75] and for elaborating different methods of quantization 
for gravity [581 M, 121 EQl EH E2] . 

3.2.2 Nonholonomic variables and fractional gravity 

Introducing the fractional canonical d-connection a D into the Einstein 
d-tensor (|22|) , following the same principle of constructing the matter source 
a T/35 as in general relativity but for fractional d-connections, we derive 
geometrically a fractional generalization of N-adapted equations ()25[) , 

"E 05 = a Ypg. (27) 

Such a system can be restricted to fractional nonholonomic configurations 
for a V if we impose a fractional analog of constraints of type (|26h 

a l% = a e a ( a N<), a C) h = 0, a n a jt = 0. (28) 

There are not theoretical or experimental evidences that for integer di- 
mensions we must impose conditions of type (128p . Nevertheless, we shall 
consider them in section [5] for deriving fractional black hole solutions which 
will mimic maximally similar ones in general relativity. 

3.2.3 Separation of equations for fractional and integer dimen- 
sions 

One of the main purposes of this work is to prove that a very general 
ansatz of type (fT6|) defines exact solutions of the fractional Einstein equa- 
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tions (|27p metric. Let us consider a fractional metric 





a Vi(x k ,v) °g l (x k ,t) a dx i ® a dx l 
+ a Va (x k ,v) °h a (x k ,v) a e a ® a e a , 
a dv + a vf(x k ,v) °Wi(x k ,v) a dx l , 
a dy A + a r]t(x k ,v) %ni(x k ,v) a dx' 



(29) 



when the coefficients will be defined below. For simplicity, we shall work 
with the "prime" dimension splitting of type 2 + 2 when coordinated are 
labeled in the form u 13 = (x^,y 3 = v,y A ), for = 1,2. This ansatz has 

one Killing symmetry because the coefficients do not depend explicitly on 
variable y 4 . 

In brief, we can write such the metric ()29f) in the form 



g = a gij a dx i ® a dx j + a h ab ( a dy a + a N% a dx k )®( a dy b + a N b k a dx k ), 



where a gij = diag[ a gi = a rji %9i] and a h ab = diag[ a h a = a rj a %h a ] and 
a N% = a Wi = a rjf %Wi and a N% = a ni = a rjf %n,i. The gravitational 
'polarizations' a 7] a and a r]f determine fractional nonholonomic deforma- 



tions of metrics, «g=[ % 9i , %h a , %N%] -> «g=[ %g % , %K, °iV fc °]@_ The 
solutions of equations will be constructed for a general source of typa 10 ! 

a T Q /3 = dm 5 [ a T 7 ; a T 1= a T 2 = a T 2 {x k ,v); a T 3 = a T 4 = a T 4 {x k )} 



Einstein ()22jl d-tensors corresponding to ansatz (|30j) reduces the Einstein 

9 Such transforms of geometric objects (with deformations of the frame, metric, connec- 
tions and other fundamental geometric structures) are more general than those considered 
for the Cartan's moving frame method, when the geometric objects are re-defined equiv- 
alently with respect to necessary systems of reference. 

10 such parametrizations of energy-momentum tensors are quite general ones for various 
types of matter sources, which (in this work) are generalized for fractional distributions; 
they can be defined by corresponding frame transform 

11 we omit such a cumbersome calculus which is similar to those presented in Refs. 
[571 1551 l6ll 170] ; for the formulas considered in this work, we have to change usual partial 
N-adapted derivatives into fractional ones considering that transform of type 1)180 might 
be performed in order to take into account the advantage that the action of Caputo 
derivative is zero for some constant coefficients 



(30) 
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equations (j27[) to this system of partial differential equations: 

1 



ar>2 
K 2 



2 a 9l a g 2 



92 



2 a g 2 

a f54 

Ha — 



9lf , an a 9'i a 92 
9i 



a 9l a 9l 

2 a gi 



'5i 



2 « 52 2 



1 



2 Q /i 3 a /t 4 1 



( a hlY a h* 3 a h* 



2 a h 4 



a '"4 1 _ 

2«/i 3 J ~ 2 ' 



1 4 



(32) 



(33) 



2 a /l 4 



2 Q /l 4 



2 a h* 



(34) 



+ 



/ ^ a h 3 d k a h A 



a R 



4 a h A 
*h 4 



Ik 



2 a h* 



a ** 
n k 



+ 



II, q 

% ^"2 ^ 4 



fc3 *; a h* 



2 Q /i 4 



0, 



2 a h* 



0, 



(35) 



In brief, we wrote the partial derivatives in the form 

a a a a a a 

a a' = 9 x a = lX i d x i a a, a a' = d 2 a = lX 2 d x 2 a a, a a* = d v a = lV d_ v a a, 

see the left Caputo fractional derivatives ([8]). 
Configurations with fractional Levi-Civita connection a V can be ex- 
tracted by imposing additional constraints 

e^ln A14 , efc Wi — e; w k , 



(36) 



*n* = 0, d { a n k = d k a rii 

satisfying the conditions ([281) . 

Following the method considered in [70] for integer dimensions, we can 
construct 'non-Killing' solutions depending on all coordinates when 

Q g = a 9i {x k ) "dx* ® a dx i + a io 2 (x j ,v,y A ) a h a (x k ,v) a e a ® a e a , 



a„3 



a dy 3 + a Wi {x\v) a dx\ a e 4 = a dy A + a n i (x k , v) a dx\ (37) 



for any a oj for which 

a e fc a uj = d k a oj+ a w k a u* + a n k d y 4 a u = 0, 

when (|37p with Q w 2 = 1 is of type (|30p . The length of this paper does not 
allow us to study such general fractional solutions. 
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4 General Solutions for Fractional Einstein Equa- 
tions 

There is a very important separation property in the system (|32|) — f)35|) 
which allows us to find exact very general solutions for such equations both 
for non-integer and integer dimensions. For instance, if the coefficient 
a gi(x k ) is known, we can find from (132p the value of a g2{x k ) (we may 
consider an inverse situation when a gi(x k ) is to be computed for a known 
value of a g2(x k )). Similarly, we can define from (|33l) the value of a h 3 (x k ,v) 
from a h 4 (x k ,v), or inversely. Having defined a h a , at the third step, we 
compute the N-connection coefficients: a Wi(x ,v) are certain solutions of 
algebraic equations (f34"|) . Then, we compute a rii(x k ,v) by integrating two 
times on v in (j35j) . Such a property of exact integration of the field equations 
can not obtained for fractional gravity models with the RL derivative but 
exists for those based on the Caputo one, correspondingly adapted to some 
classes of nonholonomic distributions. 

The explicit form of solutions of the fractional Einstein equations depend 
on the values of the coefficients in the ansatz for metric and source. Let us 
show how such solutions are constructed for ansatz of type (|30p . 

4.1 Solutions with a h* 3>i ± and a T 2A ^ 
Such metrics are defined by ansatz 

a g = e a ^ xk) a dx i ® a dx i + h 3 {x k , v) a e 3 ® a e 3 + h 4 (x k ,v) X® X, 
Q e 3 = a dv+ a Wi (x k ,v) a dx\ a e 4 = a dy A + a m{x k ,v) a dx i (38) 

with the coefficients being solutions of the system 

= 2 a T 4 {x k ), (39) 
a h\ = 2 a h 3 a h 4 a T 2 (x\v)/ a (j ) *, (40) 

a /3 a wt + a a t = 0, (41) 
X* + °7 X = 0, (42) 

where a (j) = ml , — I, a j = (\n\ a h 4 \ 3/2 /\ a h 3 \Y , (43) 

V\ Q/l 3 a h 4 \ ^ > 

a ai = a h\d k a cf>, a (3 = a h\ a 4>* . 
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For a h\ / 0; a T 2 + 0, we have a (p* / 0. The exponential function e 

in (j38|) is the fractional analog of the "integer" exponential functions and 

called the Mittag-Leffer function E a [{x— x) a ]. For a tp(x) = E a [{x— 1 x) a ], 

a 

we have d_iE a = E a , see (for instance) |42| . For simplicity, hereafter we 
shall write usual symbols for functions as in the case of integer calculus, 
but providing a label a considering such fractional construction as certain 
Taylor series in [18j. 

It is possible to consider any nonconstant a (p = a 4>(x l ,v) as a generating 
function, we can construct exact solutions of (f39]) - (j4"2|) . We have to solve 
respectively the two dimensional fractional Laplace equation, for a gi = 
°92 = e • Then we integrate on v, in order to determine a h^, "h^ 

and a rii, and solving algebraic equations, for a Wi. We obtain (computing 
consequently for a chosen a (f)(x k ,v)) 



'91 = 


92 — e. , IT-3 


7i 4 = 


OCT ( k\ _i_ ri r 






Wi = 


-1 <W V, 


"rtj = 





a <j)*{x\v)\ 
Q (/.(x fc ,u)])* 



(44) 



where o^4(x fc ), "n^ (x*) and 2 n k ( x% ) are integration functions, and 

a 

lV I v is the fractional integral on variables v. 

Here we note that the solutions (j44p . and, in general, almost all solutions 
in fractional calculus with left Caputo derivatives, can be considered as some 
series expansions and relevant fractional differential equations as is sketched 
in ref. [18]. In such an approach, for various types of fractional functions, 
there are defined a kind of fractional Taylor series of infinitely fractionally- 
differentiable functions. 

To construct exact solutions for the Levi-Civita connection "V, we have 
to constrain the coefficients (|44p to satisfy the conditions (|36p . For instance, 
we can fix a nonholonomic distribution when 2 n k (%) = and \n\~ {xf) 

a a 

are any functions satisfying the conditions 3j fn^ (xj) = dj. frii [xA . The 
constraints on a (p(x k ,v) are related to the N-connection coefficients a Wi = 

a 

— 9j a (f)/ a (f)* following relations 

( a Wi [ a </>])* + a Wi[ ( °M a <P\T+h a h 4 [ a cf>] = o, 

d i a w k [ a <t>] = d k a Wi [ a 0], (45) 
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where, for instance, we denoted by Q /i4[ a 4>] the functional dependence on 
a 4>. Such conditions are always satisfied for a <f> = a <f)(v) or if a cj) = const 
when a Wi(x k ,v) can be any functions as follows from (|4ip with zero a (3 
and ° l aii, see (|4"3"]l ). 



4.2 Three other important classes of solutions 

If any of the conditions a h^ 4 ^ is not satisfied, we can construct 
another types of solutions for certain special parametrization of coefficients 
for ansatz (|44p subjected to the condition to be a solution of equations 

(ESHM]). 



4.2.1 Solutions with a h\ = 

The equation ([33]) can be solved for such a case, a h\ = 0, only if a Y2 = 
0. Any set of functions a Wi(x k ,v) can define a solution of (p4l) . and its 
equivalent ([H]). because the coefficients Q /3 and Q aj, see (j4*3"p . are zero. 
The coefficients Q nj are determined from (I42p with a h\ = and any given 

OS 

Q /i3 which results in a n k = "n^ + 2 n fc a h%- It is possible 

to take a g\ = a g2 = e "iK 35 ) ; with a ip{x k ) determined by (I39p for a given 
Q T 4 (x fc ). 

This class of solutions is given by ansatz 

°g = e "cte' <g> <W (46) 

+ Q ^ 3 (^ fc , v) a e 3 <g> a e 3 + ^/i 4 (x fc ) a e A ® Q e 4 , 
Q e 3 = + a Wi(x k ,v) a dx\ 

1n k (x l ) + %n k (x l ) lV I v h 3 a dx l 



Q e 4 = a dy i + 



for arbitrary generating fractional functions a h^(x k ,v), a Wi(x k ,v), Qh4(x k ) 
and integration fractional functions (x 4 ) and 2 n k {%) ■ 

A subclass of solutions for the Levi-Civita connection can be selected 
from (I46p by imposing the conditions (136|) 



2«fc (a?*) = and d { "n k = d k 
a w* + di % = and % a w k = d_ k a w u 
for any such a Wi(x k , v) and Qh4(x k ). 
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4.2.2 Solutions with a h* 3 = and a h\ + 

The ansatz for metric is of type 

» g = e a ^ k ) a dx i <g> a dx i 

- %h 3 (x k ) Q e 3 Q e 3 + a h 4 (x k ,v) a e% Q e 4 , 
Q e 3 = a dv+ a Wi (x k ,v) a dx\ (47) 
Q e 4 = a dy 4 + a m{x k ,v) a dx\ 

where a g\ = a g2 = e a ' t t'( xk ) j for a ip(x k ) being a solution of (j39j) with a 
given a Ti{x k ). A function a hi(x k ,v) solves the equation PU|) for a h 3 = 0, 
which can be represented 

° h T ~ ~ 2 ohs a h 4 a T 2 (x k ,v) = 0. 

The solutions for the N-connection coefficients are 

"Wi = -1; a 4>/ a 4>*, 

a n t = ?n*(x*) lB /«[l/(vT^) 3 ], 

when Q ^ = ln| a /i 4 /VI o^3 "Ml- 

The Levi-Civita conditions (|36p for (|47|) are 

a a 

1 n fc (a?*) = and ^ = d fc " n i; 
( + a w t [ a 0] ( °7i 4 [ Q ^])* +| = 0, 

For small fractional deformations, it is not obligatory to impose such con- 
ditions. We can consider integer Levi-Civita configurations and then to 
transform them nonholonomically into certain d-connection ones. 

4.2.3 Solutions with a 4> = const 

Fixing in (|43p a (j) = a (fro = const and considering a h 3 7^ and a h\ 7^ 
0, we get that the general solutions of (f39]) - (j4"2|) are 

a g = e « dx i q a dx i _ 

%h 2 { a f*{x\v)} 2 \ ( x \v) I a e 3 <g> a e 3 + a f 2 (x\v) Q e 4 ® Q e 4 , 
°e 3 = a dv+ a Wi {x k ,v) a dx\ a e 4 = a dy 4 + \(x>) a dx\ (48) 
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where q h = const and a g\ = a g2 = e "'M 3 ' ) ; with a ip(x k ) being a solution 
of (|39p for any given a T4(x k ). 
Using the fractional function 

(x» = % [0] (*<) - lT 2 (x fc ,^)[ «/ 2 (*»] 2 , 

we write the fractional solutions for N-connection coefficients a N? = a Wi 
and a Nf = a n; in the form 



Wi 



d i a ^ r (x k ,v)/ a ^(x k ,v) (49) 



and a n k = ?n k {x l ) + %n k (s*) % I " f * ^ ,V }} 2 a ?T (z>) .(50) 

If (a: 4 , u) = ±1 for a T 2 — > 0, we take a ?4[o] (x l ) = ±1. For such con- 
ditions, the functions a /i 3 = — q/i 2 [ a f* [x\ v)] 2 and a /i 4 = a f 2 (x % ,v) 
satisfy the equation ([30]) . when ^/| a /i 3 | = °/i(y / | a /i4|)* is compatible with 
the condition a </> = Q </>o. 

The subclass of solutions for the Levi-Civita connection with ansatz of 
type ([38]) is subjected additionally to the conditions ([36]) . in this case on 
fractional function a ?x- For instance, we can chose that 2 n k ( x% ) = and 

a a 

fn k (x l ) are any functions satisfying the conditions <9j fn k = d_ k \n%. 

a 

The constraints on values a Wi = — °?x/ °<?x resu lt in constraints on a ?x, 
which is determined by Q T2 and a f, 



a 



( a w t \ a ^\y + a Wl \ \ r ] ( a h A \ \ T \r +d i a h 4 [ a , r ] = o, 

di a M^r]=d k a Wi[ a ^}, (51) 

where, for instance, we denoted by a h^ a ?x] the functional dependence on 
a Sx- Such conditions are always satisfied for cosmological solutions with 
a f = a f( v )- For "D, if a T2 = and a (j) = const, the coefficients 
a Wi(x k ,v) can be arbitrary functions. The simplest parametrization are for 
Q ?x = 1; this does not impose a functional dependence of a W{ on a sx) as 
follows from (]4ip with zero a f3 and Q aj, see (I43p . To generate solutions 
for a V such a Wi must be additionally constrained following formulas ([51 p 
re-written for a Wi[ a sr] — > a u>i(x fc ,v) and a /i4[ a sx] - >• (a; 2 ,^) • 

Finally, we emphasize that any solution Q g = { a g a '/3'(u a ')} of the Ein- 

stein equations (|25p and/or (|24|) with Killing symmetry (for local co- 
ordinates in the form y 3 = v and y 4 = y) can be parametrized in a form 
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derived in this section. Using frame transforms of type a e a = e a a a e a i, 
with a g a j3 = e a ' a e ^ a g a 'i3', for any Q g a/ 3 (f29j) . we relate the class of such 
solutions, for instance, to the family of metrics of type (j38|) . 

5 Fractional Spacetimes and Black Holes 

A fractional spacetime is with very different rules of computing local par- 
es 

tial derivatives (via integration) of type ix id i ^fy instead of usual partial 
derivatives <9j. The actions of such operators are very different on singu- 
lar functions, for instance, parametrized in the form ([I]). So, the singular 
etc structure of solutions of fractional Einstein equations ([25]) and/or (f24"|) 
should differ substantially from that for integer dimensions in general rela- 
tivity. 

There is a series of very important questions to be solved in fractional 
models of gravity: 1) for instance, if such theories contain black holes so- 
lutions? 2) if such solutions for fractional black holes can be constructed, 
what are their properties? 3) what may happen with an integer dimen- 
sional black hole into a fractional background, for instance, of solitons? We 
should also analyze the problems: 4) what may happen with integer and 
non-integer dimensional black holes under fractional/nonholonomic Ricci 
flows? 5) what is the status of singular solutions and horizons in fractional 
classical and quantum gravity etc. 

In this section, we prove that black holes really exist in fractional gravity 
contrary to the hope that involving a new type of derivative calculus ((2|), and 
changing respectively the differential spacetime structure, we may eliminate 
"ambiguities" with singularities etc. The concepts of black hole, singularity 
and horizon seem to be fundamental ones for various types of holonomic and 
nonholonomic, commutative and noncommutative, pseudo-Riemanann and 
Finlser like, fractional and integer etc theories of gravity. We also provide 
solutions for the questions 1-3) above and leave 4) and 5) for our future 
investigations. 

5.1 Fractional deformations of the Schwarzschild spacetime 

We consider a diagonal integer dimensional metric e g depending on a 
small real parameter 1 > e > 0, 

£ g = -d£®d£-r 2 (£) cM<g>cM-r 2 (£)sin 2 i? dip ® dip + w 2 {£) dt® dt. (52) 
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The local coordinates and nontrivial metric coefficients are parametrized: 



x 
91 



92 



#,y =v = <p,y 



for £ 



dr 



2^o 

r 



1/2 



= t, (53) 
-r 2 (£)sin 2 i?, h 4 = w 2 {£), 

2^o 

r 



and w 2 (r) = 1 



e 

72 - 



For e = in variable £(r) and coefficients, the metric (|52|) is just the the 
Schwarzschild solution written in spacetime spherical coordinates (r, t) 
with a point mass //o- 

We search for a class of exact fractional vacuum solutions of type (|29|) 
when the fractional metrics are generated by nonholonomic deformations 



'9i. 



x f]igi and h a = a rj a h a and some nontrivial 



Wi, 



*rii, where (gi,h a 



are given by data (|53|) and parametrized by ansatz 

= - a m(t,$,0) a d£® a dt-m(ti^,0)r 2 (O a <M® a d$ (54) 



when the coefficients will be constructed determine solutions of the system 
of equations (p]) - dl2j) with a T^ = 0. 

The equation (|40p for Q T2 = is solved by any 

(55) 



*h 3 
*h 4 



a b*) 2 



a, 2 



r/ 3 (£,tf,<^)r 2 (£)sin 2 tf, 
m {t,d,^0)m 2 {i), 

2 



Qh) 2 \hi/h^\ ( a r]i\\ . We consider q/i = const (it must 



q ii — 2 in order to satisfy the condition (I40p with zero source), where "774 



for 

be %h 

can be any function satisfying the condition a r/| 7^ 0. This way, it is possible 
to generate a class of solutions for any function a 6(£, (p, 9) with a b* 7^ 0. 
For classes of solutions with nontrivial sources, it is more convenient to work 
directly with fractional polarizations "774, for a r}\ 7^ 0. In another turn, 
for vacuum configurations, it is better to chose as a generating function, for 
instance, a h 4 , for a h\ 7^ 0. The fractional polarizations a r\\ and a r\i, when 
a f]i = a r]2r 2 = e ^(£>*) (for fractional configurations, this is not just the 
exponential function but its fractional version), from (|39p with a T4 = 0, 
i.e. from a ip" + a ip" = 0. 
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Putting the above coefficient in (|54p . we construct a class of exact vac- 
uum solutions in fractional gravity defining stationary fractional nonholo- 
nomic deformations on a small parameter e of the Schwarzschild metric, 



■g 



a S<p 
a 5t 



-4 



2 

,2/ 



(56) 



a V*(S,0,<P,9)™ 2 (O a $t® a 5t, 



In general, the solutions for fractional metrics (|54p . or (|56p . do not define 
black holes and do not describe obvious physical situations. We can con- 
sider that for general nonholonomic fractional deformations a usual " integer" 
black hole may " dissipate" into a fractional structure of a more sophisticate 
fractional spacetime ether with a not defined status of singularities of coef- 
ficients of metric. 

In next subsections, we show that it is possible to chose certain non- 
holonomic distributions when the singular character of the coefficient tu 2 (!;) 
vanishing on the horizon of a Schwarzschild black hole result in physical 
properties of usual black holes. 



5.2 Non— integer gravitational ellipsoid configurations 



It is possible to extract a class of metrics (J56J) defining fractional de- 
formations of the Schwarzschild solution depending on parameter e with 
possible physical interpretation of fractional gravitational vacuum configu- 
rations with spherical and/or rotoid (ellipsoid) symmetry. We chose in (I55p 
a generating functions of type 



a b 2 



(57) 



considering, for simplicity, only linear decompositions on a small parameter 
e. For ([ST]) , we get 



So, we can compute in ([56]) the coefficients a h^ and a /i4 and corresponding 
polarizations "773 and a rjA. using formulas (|55p . Here we note that if we put 
e = we can generate fractional deformations of the Schwarzschild solution 

2 r 12 

not depending on parameter a, when a b = a q and ( a b*) 
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Fractional rotoid configurations are generated for 



> = 1 - and «s = fP, S in(uw + (58) 

r 4 a [i 

when a /i(£, i?, if) = no + 92) describes fractional locally anisotropic 

polarized mass. The constants /ioi^o and (fo and arbitrary functions 

i?, V?) and go( r ) have to be defined from some boundary conditions, 
with e treated as the eccentricityr^l of an ellipsoid. 

The fractional stationary rotoid solutions for the Schwarzschild metric 
in general relativity can be written in the form 

-4 [(v^i)*] 2 [i + e( a Vv^i)*/(\^)*] a <ty>® a ty 

+ ( a q + e a s) a 5t® a 5t, (59) 
a 5(p = a d<p + a Wl a d£ + a w 2 a dti, 
a 5t = a dt + f m Q d£ + ?n 2 a di?, 

where functions a g(£, 1?, </j) and a s(£, ip) are given by formulas (|58|) . The 
N-connection coefficients, a Wj(£, 1?, <£>) and a n« = are subjected 

to conditions of type ([36]) . 



a w 2 (In I Q «; 2 |)* = a wl - a w[, a w* / 0; 

or a w' 2 - a w[ = 0, a w* = 0; 
?ni(£>0)- ?n5(£,0) = 

and a tp(^,'&) being any function for which "-i/;" + "-(//' = 0. 

For small eccentricities, a metric (|59|) defines stationary fractional con- 
figurations for so-called black ellipsoid type solutions. Their stability and 
properties can be stated and analyzed by adapting to fractional calculus the 
methods elaborated in [52J, [531 G3] ( a summary of results and generalizations 
for various types of locally anisotropic gravity models in Ref. [55J, see simi- 
lar constructions for a noncommutative Finsler version of rotoid spacetimes 
in [73]). 



12 e is an eccentricity because, for integer configurations, the coefficient hi = b 2 = 
V 5 ) ^(f) becomes zero for data ((58} if r+ ~ 2/io/[l + s ^-r sin(uoy + yo)]; 
this is the "parametric" equation for an ellipse r+Up) for any fixed values , uo, yo and 
Mo 
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5.3 Stationary black holes in fractional solitonic configura- 
tions 



Finally, we analyze two types of fractional rotoid solutions (in particular 
the Schwarzschild one) imbedded into 1) integer solitonic background and 
2) into a fractional solitonic background denned by an exact solution of a 
fractional differential equation. 



5.3.1 Integer solitonic background 

There are static three dimensional solitonic distributions fy(£, i?, <p, 8), 
defined as solutions of an (integer) solitonic equation 

77" + e(rj' + 6i] 7]* + rj***)* = 0, e = ±1, 

resulting in stationary black ellipsoid-solitonic fractional metrics generated 
as further deformations of a metric fot £ 



g ([59]) 1^1 The generated fractional 



solitonic gravitational metrics are of type 



x 1p ( ot 



l + e- 



*cW<8> a d$) 
1 / 



a 5p 
a 5t 



+i 1 ( a q + E a s) a 5t ® a 6t, 
a dip + a Wl a d£ + a w 2 a d$, 
a dt + f m a di + > 2 a d$, 



(60) 



where the N-connection coefficients are taken the same as for (|59p. 

The metrics (j60j) are of type (J56j) . So, they positively define vacuum 
solutions of fractional Einstein equations. 



5.3.2 Fractional solitonic backgrounds 

An interesting property of the solutions (j!5|) is that they depend on a 
generating fractional function a f 2 (x l ,v) which is a general one, but con- 
strained to some conditions in order to generate a corresponding class of 
solutions. For instance, we can consider that the fractional metric (|60p is 
additionally deformed nonholonomically by a solution a p(v), when a p* = 

13 a function rj can be a solution of any three dimensional solitonic and/ or other non- 
linear wave equations 
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d v p = it) d v a p and iv d v being the left RL derivative on v, of a fractional 
differential equation 



M V) + ^(v) v 



2 z(v) 



(61) 



for some suitable functions 1 z{v) and 2 z{v). This fractional equation can 
be solved in general form [BJ, 



*p(v) = 1V I V 

p=0 



2 z(a) 



+ ^(v- V) + 2 c, (62) 



where x c and 2 c are constants. If 1 z(v) = and 2 z(v ) = A a p(v), where A € 
R, the equation (|6ip transform into a fractional Euler-Lagrange equation 
for the one-dimensional fractional oscillator. Its solution is different from 
(|62p but also can be expressed as a series (see details in [6]). 

We construct a fractional black rotoid solution embedded both into an 
integer solitonic background and a one-dimensional fractional gravitational 
generalized oscillator (|62j) if we consider the (vacuum solution) metric 



g 



-e ( a d£ ® a d£ + a cM a d$) - 4 [(VI?? a tf a p\T 



l + e- 



a 5<p 
a St 



(V\v a Q a p\Y \V\v a Q a p\ 

+ V a p( a q + e a s) a 6t® a 5t, 
a dp + a wi a d£ + a w 2 a d$, 
a dt + f ni + in 2 a d$, 



where the N-connection coefficients are taken the same as for (|59p . This type 
of metrics are also of type (|56p , with a different nonholonomic structure, and 
also define vacuum solutions of fractional Einstein equations. 

The fractional nonholonomic black hole type solutions (|59|) and (|6U|) , and 
related solitonic and fractional oscillator extensions, are stationary. 



6 Discussion, Conclusions and Further Develop- 
ments 

Different concepts of dimension and related definitions of derivatives 
(Riemann-Lesbegue, Caputo etc) and integrals (Lesbegue-Stieltjes, Riemann- 
Stieltjes, fractional integrals etc) were introduced long time ago and studied 
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at present in various branches of science. In this work, we developed a geo- 
metric approach to fractional gravity theories based on fractional nonholo- 
nomic manifolds, nonlinear connections and generalized solutions of frac- 
tional Einstein equations, deriving all constructions from the Caputo left 
fractional derivative. Following such a direction, we preserve a maximally- 
possible analogy between spaces of integer and non-integer dimensions. For 
study of general dynamical and evolution theories of gravitational and mat- 
ter fields (classical and quantum models), we have to apply various meth- 
ods from the geometry of nonholonomic manifolds (originally elaborated for 
Lagrange-Finsler spaces) developed for a correspondingly adapted fractional 
calculus. 

In our approach, the fractional gravity can be naturally defined on frac- 
tional manifolds following the same principles as in general relativity. More 
than that, there is an unified method allowing to derive both types of metric 
compatible generalizations of the Einstein theories, with a fractional version 
of the Levi-Civita connection, and of generalized Lagrange-Finsler models, 
with the canonical/Cartan distinguished connections. The surprising result 
is that the fractional Einstein equations can be integrated in general form 
by adapting the constructions with respect to certain classes of frame elon- 
gated linearly by a nonlinear connection (N-connection) structure. Both 
for the fractional and integer dimension (pseudo) Riemannian spaces, the 
coefficients N-connections are defined by certain off-diagonal terms of met- 
rics. For Lagrange-Finsler theories, the N-connections are determined as a 
fundamental geometric structure induced by a fractional/integer generating 
function. 

In this work, we classified possible integral varieties of fractional nonholo- 
nomic gravitational field equations. We also provided a study of Schwarz- 
schild type black holes and rotoid configurations in fractional spacetimes. 
The main conclusion was that the black holes exist also in fractional gravity 
and that such configurations may survive for certain types of small non- 
holonomic fractional deformations of metrics and in integer, or fractional, 
solitonic backgrounds. Nevertheless, more general classes of fractional non- 
holonomic transforms result in exact solutions with not obvious physical 
interpretation. 

There are many directions of investigation left to explore in fractional 
calculus and geometry and applications in various sciences. Here we outline 
seven important ones being related to possible developments of our former 
results in (integer dimensional) geometry of nonholonomic manifolds/ bun- 
dles, and applications to modern classical and quantum physics theories: 
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1. Foundations of fractional classical and quantum field theories with 
possible noncommutative, supersymmetric, string, brane generaliza- 
tions etc. There is a number of publications in such directions [H 
El [81 O HH [HI [13 [231 SD S21 ETl [89]. Our proposal is to analyze if 
such constructions can be redefined for nonholonomic fractional non- 
commutative/ supersymmetric etc spaces and string/brane models, 
see reviews of such results in our works [46 \ l4"71 l4"9l 1821 [55] , when for 
adapted nonholonomic distributions the constructions the fundamen- 
tal field equations became integrable in very general forms. 

2. A theory of fractional nonholonomic Ricci flows was proposed recently 
in [75] . This generalize in a fractional fashion a series of our papers 
on nonholonomic Ricci flows of Einstein, Lagrange-Finsler and non- 
commutative/ nonsymmetric spaces |62l l62| [68] and exact solutions 
of the Einstein equations and Hamilton's equations and applications 
in cosmology and astrophysics [651 EU [67] . Solutions with fractional 
evolution of geometric objects present a strong theoretical arguments 
for study such theories. 

3. Exact solutions in fractional gravity and applications in modern cos- 
mology and astrophysics. The first examples of such fractional exact 
solutions were constructed and analyzed in previous two sections of 
this paper. We proved that our version of fractional Einstein equa- 
tions can be solved in general form and analyzed models of black 
hole/ellipsoid fractional metrics. Here we note that there are various 
attempts to relate modern cosmological experimental data to effects of 
fractional/fractal gravity and matter field interactions [27 ] l2"8 [ [35 ] l8[ 19] . 
With respect to our former works on general/exact solutions in various 
models of gravity theories and applications, we suppose to be of inter- 
est in modern mathematical physics and gravity such issues: a) Exact 
solutions for models of fractional gravity of arbitrary dimensions and 
with different types of connections generalizing [TQl [551 [57] ; b) locally 
anisotropic fractional cosmology models |74[ 155]: c) various fractional 
brane, noncommutative, anisotropic Taub NUT, wormhole and black 
hole solutions [731 [801 Ell [52l [531 [Ml ESI E3] . 

4. Fractional Clifford spinor structures and Dirac operators have been 
considered, for instance, in |32[ 124]. Following an "exactly integrable" 
fractional nonholonomic approach, we find such perspective for fur- 
ther developments on fractional Clifford geometry and applications. 
There are at least five different sub-directions: a) fractional Clifford- 
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Lagrange/-Finsler/ - Hamilton spaces and higher order generaliza- 
tions, see original results for integer dimensional spaces in [44\ [4~8"1 
[MIES]; b) fractional Clifford-Finsler algebroids, extending [56J, and 
c) nonholonomic non-integer gerbes as generalizations from [78]; d) 
exact solutions of fractional Einstein-Dirac equations adapted to N- 
connection structure as we considered in [79 [ [83] and e) fractional 
models of noncommutative geometry, Dirac operators and noncommu- 
tative gravity and field interactions, see original results in [681 EH EI] 
and Part III of [55]. 

5. Fractional Lagrange-Finsler and Hamilton- Cartan spaces and higher 
order generalizations. This direction is originally considered in [1] as 
a fractional generalization (using the Riemann-Liouville integral) of 
some results from |2H 120]. Here we cite our alternative/ complimen- 
tary constructions with higher order (super) vector bundles/ Clifford 
bundles and super-strings [Ml EH SHI E21 EB H7J H2] which can be gen- 
eralized for fractional higher order calculus, adapted to nonholonomic 
distributions and fractional Caputo derivatives, with exactly integrable 
string/brane/ spinor and modified gravitational field equations. 

6. Fractional quantum gravity, deformation quantization, renormalizabil- 
ity, and gauge models. Such theories, for instance, in fractional ver- 
sions, are supposed to provide new tools in quantization of gravity and 
matter field interactions, see [U [HI [T71 [381 ES] . To elaborate a geomet- 
ric approach to fractional quantum theories, including quantum grav- 
ity, is a topic for future investigations. Working with nonholonomic 
distributions, such constructions can be derived fractionally from our 
former results on deformation quantization of gravity and Lagrange- 
Finsler and Hamilton-Cartan spaces |58[ IBTil [2], A-brane fractional 
quantization [M] and bi-connection quantization [7TJ [72] 

7. Diffusion, kinetics thermodynamics and fractional nonlinear dynamics 
with solitons, chaos and fractals. This is, perhaps, the most elabo- 
rated direction in fractional calculus with a number of applications, 
see QUI [861 dH EH E3 [88] and references therein. Such fractional 
constructions adapted to nonholonomic distributions are possible for 
locally anisotropic (in general, supersymmetric) stochastic calculus, 
see [45] and Chapters 5 and 10 in monograph [49|, and kinetic theories 
and anisotropic thermodynamics [501 [5T] . There is also recent our pa- 
pers with bi-Hamilton structures, solitonic hierarchies and encoding 
the solutions of gravitational field equations in Einstein-Finsler gravity 
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and Ricci-Lagrange theories [3j [691 EH] ■ Solitonic hierarchies encoding 
of solutions of fractional systems of equations related to fundamental 
theories present a substantial interest for fractional dynamics, evolu- 
tion and field interactions. 

A series of our further works will be devoted to researches concerning 
above directions. 
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